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The coupling of a young stellar disc with the molecular 
torus in the Galactic centre 
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ABSTRACT 

The Galactic Centre hosts, according to observations, a number of early-type stars. 
About one half of those which are orbiting the central supermassive black hole on 
orbits with projected radii > 0.03 pc form a coherently rotating disc. Observations 
further reveal a massive gaseous torus and a significant population of late- type stars. In 
this paper, we investigate, by means of numerical iV-body computations, the orbital 
evolution of the stellar disc, which we consider to be initially thin. We include the 
gravitational influence of both the torus and the late-type stars, as well as the self- 
gravity of the disc. Our results show that, for a significant set of system parameters, the 
evolution of the disc leads, within the lifetime of the early-type stars, to a configuration 
compatible with the observations. In particular, the disc naturally reaches a specific - 
perpendicular - orientation with respect to the torus, which is indeed the configuration 
observed in the Galactic Centre. We, therefore, suggest that all the early-type stars 
may have been born within a single gaseous disc. 
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Over the past two decades, nearly 200 early-type stars 
have been revealed in the innermost parsec of our Galaxy 
(see Genzel et al. 2010 for the most recent review ; 
Allen et al.lll990l:lGenzel et al.ll2003) : | Ghez et al|2003l . 120051 ; 

20091 . |201oh . Observa- 



Paumard et al 



2006; 



Bartko et al 



tions suggest that these stars are orbiting a highly con- 
centrated mass, which is associated with the compact ra- 
dio source Sgr A*. It is widely accepted that this source is 
powered by a supermassive black hole (SMBH). Its mass 
and distance from the Sun are estimate d to be approxi 



mately 4 x 10 6 Mg and 8 kpc, respectively dGhez et al. 2003; 
Eisenhauer et al.ll2005l ; iGillessen et al1l2009al lbl lYelda et all 



20101 ). 



Accordin g to the most recent observations of 
iBartko et ail l|2009l , l2010t ), the majority (136) of the 
early-type stars observed in the Sgr A* region are located 
at projected distance 0.03 pc < r < 0.5 pc from the SMBH. 
Roughly one half of these stars appear to form a coherently 
rotating disc- like struct ure, the so-called clockwise system 
(CWS; discovered by iLevin fc Beloborodovl l2003h . The 
remaining stars are randomly scattered off the CWS plane. 
Nevertheless, some authors report the existence of a second 
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coherent structure at a large angle with respect to the CWS 
- the counte r clock wise system (CCWS; first mentioned by 
IGenzel et all l2003h . Given these two discs, a significant 
number of the early-type stars are still not belonging to 
either of the structures. 

Observations have further established that all the early- 
type stars between 0.03 pc and 0.5 pc from the SMBH are 
either Wolf-Rayett stars or O- or early B-stars (WR/OB 
stars) l|Bartko et al.ll2009l . l201ol ). Evolutionary phases of in- 
dividual stars indicate that all of them have been formed 
6±2 Myr ago w ithin a short period of time, probably not ex- 
ceeding 2 Myr (Paumard ct al. 2006]). The presence of such 
stars so close to a SMBH is rather suprising. In particular, 
the tidal field of the SMBH is strong enough to prevent stan- 
dard star formation mechanisms. Hence, various hypotheses 
have been suggested to explain the origin and configuration 
of the WR/OB stars observed in the Sgr A* region. 

In situ fragmentation of a self-gravitating gaseous disc is 
probably the currently most wid ely accepted formation sce- 
nario for the stars of the CWS l|Levin fc Beloborodovl 120031 ; 
iPaumard et al.ll2006l ). This process was theoretically pre- 
dicted to form stars in a ctive galactic nuclei around SMBHs 
of masses 10 6 -10 10 M Q (ICollin fe Zahnlll999T ). However, as it 
naturally forms stars in a single disc-like structure, it fails 
to explain the origin of the stars observed outside the CWS. 
Many authors have, therefore, been seeking a mechanism 



2 J. Haas, L. Subr, and P. Kroupa 



that could have scattered these outliers from the parent disc 
plane. 



It has been shown by ICuadra et al.l ((2008) that two- 
body relaxation of the parent disc does not yield the ob- 
served large inclinations of the outliers with respe ct to the 
disc plane. According to iKocsis fc Trem ainc (2010), some of 
them may have been brought to their positions by vector 
resonant relaxation between the disc and the cluster of late- 
type st ars, which has also been reported in the Sgr A* 
regio n (|Genzel et all 120031 ; ISchodel et all 120071 ; |Po et all 
120091 ). However, it is still unclear whether this process can 
explain the origin of the stars with line-of-sight angular mo- 
menta opposite to that of the stars within the C WS. In order 
to overcome this issue, lLockmann et al.l (|2008l ) have consid- 
ered mutual interaction of two self-gravitating discs at large 
angles relative to each other. Although this mechanism in- 
deed yields the observed configuration of the WR/OB stars, 
it needs rather special initial conditions. In particular, the 
two discs must have been formed at specific angles with 
respect to each other in order to stand for the CWS and 
CCWS. Moreover, as all the WR/OB stars seem to be coeval 
(|Paumard et al.ll2006l ). the two discs must have had started 
to form stars at the same time. Since this is not very likely, 
the need of such special initial conditions represents the 
major drawback of this scenario. Simil ar scenarios, such as 
the i nteraction of two gaseous streams |Hobbs fc Navakshinl 
2009), suffer from the same problem. 



ISubr et all l|2009h and lSubrj < |2010T > have suggested that 
all the WR/OB stars in the Sgr A* region may have been 
born in a single gaseous disc. They argue that the stars 
observed outside the CWS represent the outer parts of the 
parent disc, that have been partially disrupted by the grav- 
ity of the circumnuclear disc (CND). The CND is a clumpy 
molecular tor us, that is located betwe en 1.6 pc and 2.0 pc 
from Sgr A* ^Christopher et alJfeOQSl ). The total mass of 
this structure, which is almost perpendicular to the CWS 
jPaumard et al.l 120061 ') . is estimated to be of the order of 



10 6 Mm. ISubr et all j2009h claim that the gravity of the 



CND would cause differential precession of the individual 
orbits in the parent stellar disc. Such a process would force 
the stars from the outer parts of the disc to leave the disc 
plane while the inner parts of the disc would remain un- 
touched. This core would be identified as the CWS today. 

In this paper , we further investigate the hypothesis of 
ISubr et"al] l|2009l ). In particular, we include the self-gravity 
of the parent stellar disc, and follow its orbital evolution in a 
predefined external potential by means of numerical iV-body 
computations. In addition to the SMBH and the CND, the 
external potential includes the gravity of the cluster of late- 
type stars. Even though its density profile is still unclear, 
its potential may be considered, in the first approximation, 
to be spherically symmetric, and centred on the SMBH. 

In the following section, we briefly introduce the prob- 
lem of stellar dynamics in a perturbed Keplerian potential 
and describe our model. The results of our calculations are 
presented in Section [3] and discussed in Section [4] We con- 
clude our work in Section [5] 



2 ORBITAL EVOLUTION OF A STELLAR 
DISC IN AN EXTERNAL POTENTIAL 

The gravitational potential in the vicinity of the CWS in- 
duced by the SMBH is, to a very high accuracy, Keplerian. 
Hence, it is useful to describe the individual stellar orbits 
by means of Keplerian orbital elements: semi-major axis a, 
eccentricity e, inclination i, longitude of the ascending node 
Q, and argument of pericentre u. The orbital evolution of 
the disc can then be investigated by following the evolution 
of these elements. If the only component of the external po- 
tential were Keplerian gravity of the SMBH, and if the stars 
in the disc were treated as test particles, the orbital elements 
would remain constant in time. On the other hand, with any 
additional potentials included, some of the elements may un- 
der go complex s ecular evolution. 

ISubr et alj {2009) have investigated the influence of the 
CND upon the orbital evolution of the disc. They have con- 
sidered the disc to be surrounded by the cluster of late-type 
stars, and the stars in the disc to be test particles. They 
claim that the evolution of individual stellar orbits in the 
disc is dominated by two processes: (i) precession of the as- 
cending node, and (ii) periodic oscillations of eccentricity 
and inclination (Kozai oscillations; independently theoreti- 
cally predicted by Kozai 1962 and Lidov 1962). These oscil- 
lations are diminishing with increasing mass of the cluster of 
late-type stars. They become fully negligible when the mass, 
M c , of the cluster within the radius -Rcnd of the CND fulfils 
the approximate condition M c > O.IMcnd, where Mcnd 
stands for the mass of the CND. In that case, the first time 
derivative of £1 becomes constant and can be writen as 



with 



dF 



3 cos i 1 + |e 2 



4 T K yT^ 5 



M. Rf 



Mcnd ay/GM.t 



(X) 



(2) 



where G denotes the gravitational constant and Mm rep- 
resents the mass of the SMBH. According to this formula, 
the rate of precession strongly depends upon the semi-major 
axis of the orbit. Hence, the outer parts of the disc are more 
affected by the precession than the inner parts and, there- 
fore, the disc becomes warped or, eventually, completely dis- 
rupted. 

In formula |T} and further on, we assume all angles to be 
measured in the frame where the CND lies in the xy-p\axie. 
Consequently, inclination of the orbit i is the angle between 
the symmetry axis of the CND and angular momentum of 
the star. 

With the gravity of the stars in the disc included, the 
orbital elements of individual orbits undergo random varia- 
tions due to two-body relaxation. As a result, the mean value 
of eccentricity and inclination of the orbits increases. These 
changes are, by themselves, too small to have a significan t 
impact on the overall shape of the disc (|Cuadra et al.ll2008l ) . 
However, according to formula ([T}, they accelerate the dif- 
ferential precession of the orbits. Let us, therefore, further 
investigate the differential precession of the orbits within 
a self-gravitating disc. For this purpose, we introduce the 
model of the Galactic Centre in the following way. 




Figure 1. Results for the canonical model (see Table[T]for the corresponding parameters). Only properties of the stars with m ^ 12 Mq 
are displayed. The dotted lines in the middle panels denote standard deviation for the set of 120 included realisations. Top: Evolution 
of the mean value of (left) and i (right) within different parts of the disc for one of the realisations of this model. The solid line 
describes the group of the innermost stars, followed by the dashed, dot-dashed and dotted line, which correspond to successive outer 
groups. Middle-left: Number of stars within the CWS (i. e. with angular momentum deviating from the mean angular momentum of the 
CWS by less than 20°). Middle-right: Inclination of the CWS with respect to the CND. Bottom-left: Eccentricity distribution for the 
stars after 6 Myr of orbital evolution. The empty boxes denote distribution for all the stars in the young stellar system while the grey 
ones represent only stars within the CWS. Bottom-right: Mean eccentricity of all the stars in the young stellar system (solid line) and 
within the CWS (dashed line). 
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(i) The SMBH of mass M. = 4 x 10 6 M Q is considered to 
be a source of Keplerian potential. 

(ii) The CND is modelled as a single massive particle of 
mass Mqnd orbiting the SMBH on a circular orbit of radius 
Rcnd = 1.8 pc (see Subsection 14. 21 for the discussion of this 
approximation) . 

(iii) The cluster of late-type stars is represented by a 
smooth power-law density profile, p(r) oc r - ^, and mass 
M c within the radius -Rcnd- 

(iv) The early-type stars in the disc are treated as N grav- 
itating particles, m £ [m m i n , m max ], distributed according to 
a power-law mass function dN oc m~"dm. 

The stellar orbits in the disc are constructed to be ini- 
tially geometrically circular. However, due to the additional 
spherically symmetric component of the gravitational poten- 
tial (cluster of late- type stars) , the osculating eccentricities 
do not truthfully describe real curvature of the orbits in 
space. E. g. the initial osculating eccentricity of the outer- 
most orbits in the disc is w 0.1 for M c = 0.1 and f) — 7/4. 
Initial radii of the orbits are, in accord with t he ob serva- 
tions (|Paumard et al.ll2006l ; IrBartko et al.ll2009l . |201Ch . gen- 
erated randomly between 0.04 pc and 0.4 pc. Their distri- 
bution obeys dN oc a _1 da. The disc is initially thin with 
half-opening angle Ao < 5°. The initial inclination of the 
disc plane with respect to the CND, which is defined by the 
mean angular momentum of the stars in the disc, is denoted 
Ws- 

We follow the evolution of the latter system numeri- 
cally by means of the iV-body integration code NBODY6 
l|Aarsethl [2003). The gravitational potentials of both the 
SMBH and the cluster of late-type stars have been incor- 
porated into the original code as additional external poten- 
tials. 



3 RESULTS 

The evolution of the stellar disc depends upon the shape of 
the gravitational potential, which is determined by the pa- 
rameters Mcnd, M c , P, N, m min , m max , a, A , and icws- 
Hence, in order to investigate the evolution properly, it is 
necessary to cover all the reasonable values of these param- 
eters. On the other hand, in order to demonstrate the results 
of our calculations, it is useful to define a "canonical" model 
with the parameters set to the values listed in Table [T] 

The observations sugge st that all the WR/OB stars 
have mass m > 12 M Q jpaumard et al.ll2006l ; [Bartko etajj 
l2009ll2010h . However, since it is likely that a number of un- 
detected less massive early-type stars exist in the Sgr A* 
region, we consider m £ [4Mq, 120 Mg] in the "canonical" 
model. Nevertheless, for a more convenient comparison with 
the currently available observational data, we display prop- 
erties of only a subset of the stars with mass m ^ 12 Mq in 
figures. 

Due to the stochastic nature of the studied system, the 
results should be averaged over a number of realisations with 
identical values of the model parameters in order to distin- 
guish general trends from random fluctuations. For this pur- 
pose, we first considered 120 realisations of the "canonical" 
model. It has, however, turned out that the results become 
statistically relevant already for 12 realisations. Hence, we 
consider only 12 realisations of all the other models discussed 



Table 1. Parameters of the "canonical" model (see Fig. [T] for the 
corresponding results). 



Mm = 4 x 10 6 M Q Mcnd = 0.3 M. M c = 0.03 M, 
m = 4-120 M Q a = 1 p = 7/4 

N = 200 i cws = 70° A = 2.5° 

Total mass of the young stellar disc « 6.8 X 10 3 Mq 
(> 90 % by « 140 stars with m > 12 Mq) 



in this paper in order to shorten the necessary computa- 
tional time (w 5 hours on 3 GHz CPU per run; « 2000 runs 
in total). 

Since we attempt to explain the configuration of a spe- 
cific observed system, every single realisation represents a 
possible course of its evolution. We thus show the standard 
deviation for some of the key quantities for a more thorough 
description of the set of possible evolutions. 

The results for the canonical model are shown in Fig. [I] 
The top-left panel demontrates the differential precession of 
the orbits in the disc. It shows the evolution of the mean 
value of within different groups of stars, which are de- 
termined by their initial distance from the centre. It turns 
out that the precession of the ascending node affects more 
strongly the orbits in the outer parts of the disc (dotted and 
dot-dashed lines) than those in the inner parts (dashed and 
solid lines). This result is in accord with formula (O and 
proves the gradual deformation of the disc. 

Our results further show that the precession of the as- 
cending node in the outer parts of the disc is globally accel- 
erated. We attribute this effect, which becomes significant 
on longer time scales, to the evolution of inclination due to 
two-body rel axation of t h e disc . Such an acceleration was 
not found by ISubr et "all (|2009h as they had neglected the 
gravity of the stars in the disc. 

The sudden drop of (fl) on the dotted line in the top-left 
panel of Fig.Q] is a residue of Kozai oscillations. Since the 
cluster of late-type stars is, in our canonical model, not mas- 
sive enough to suppress the oscillations of e and i entirely, 
the first time derivative of Q. also varies on the timescale of 
T K . 

The evolution of the mean inclination (i) with respect 
to the CND within different parts of the disc is shown in 
the top-right panel of Fig.[T] We can see that the inclination 
of the outer parts of the disc is decreasing (dotted and dot- 
dashed lines), while it grows and saturates at « 90° in the 
inner parts (dashed and solid lines). We further find that the 
evolution of both fi and i is similar for all the orbits in the 
inner parts of the disc. Hence, the core of the disc remains 
rather undisturbed and coherently changes its orientation 
towards perpendicular with respect to the CND. This effect 
can be seen in the top panel of Fig. [2] which shows the direc- 
tions of angular momenta of the individual stars in the disc 
after 6 Myr of orbital evolution for one of the realisations of 
the canonical model (the initial state is denoted by an empty 
circle). Our results prove that the compact group at incli- 
nation w 90° is formed by the stars from the inner parts of 
the disc, while the remaining scattered stars represent the 
entirely dismembered outer parts. Hence, we see that the 
dynamical evolution of the initially thin stellar disc leads to 
a configuration similar to that observed in the Sgr A* region 
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follow within our calculations the relative number of stars, 
Ncws/N, which belong to the CWS. As can be seen in the 
middle-left panel of Fig.[T] this number reaches, within our 
model, the value of ~ 0.5 at t = 6 Myr. 

As computed bv lPaumard et ah! l|2006l ). the normal vec- 
tors of the observed CWS and the CND fulfil news ■ "cnd = 
0.01, which corresponds to the mutual angle of 89.4°. In or- 
der to confront this feature, we investigate the evolution of 
the inclination iews of the CWS with respect to the CND. 
Our computations show that iews ~ 90° at t = 6 Myr (see 
the middle-right panel of Fig.[T]), which is in a remarkable 
agreement with the observational data. 

Finally, we investigate the eccentricity distribution n e 
within the CWS and in the whole young stellar system after 
6 Myr of orbital evolution. The corresponding histograms 
in the bottom-left panel of Fig.[T] show that in both cases 
a substantial fraction of the orbits have, in accord with the 
observations, moderate eccentricities. The mean eccentric- 
ity of the stars within the CWS (see the dotted line in the 
bottom-right panel of Fig.[TJ is then w 0.25, which is some- 
what lower than th e value 0.36 ± 0.06 recently reported by 
Bart ko~et al.l (|2009T ) . However, at the current level of accu- 
racy, the observations do not provide sufficient information 
for a reliable determination of the orbital eccentricity for a 
significant number of the WR/OB stars. Hence, the eccen- 
tricity criterion should be considered only as supplemental. 



0° 



Figure 2. Angular momenta of individual stars in the young stel- 
lar disc in sinusoidal projection after 6 Myr of orbital evolution. 
The initial state is denoted by an empty circle. Latitude on the 
plots corresponds to i while longitude is related to f2. The top 
panel shows the results for one of the realisations of the canonical 
model (only stars with m > 12 Mq displayed). For comparison, 
the bottom panel illustrates the situation with negligibly small 
mass of the stars in the disc (single mass, m = 0.004 Mq, the 
other parameters are the same as in the canonical model). 



(see lPaumard et alJl2005lBartko et"alll2009j . l20ioh . In par- 
ticular, the core of the disc can be identified with the CWS 
observed today and, at the same time, the stars from the 
dissolved outer parts can stand for the WR/OB stars found 
outside the CWS. 

In order to compare our results with the observations 
more thoroughly, we further define CWS within our model 
in the following iterative way. As the zeroth step, the CWS is 
considered to be formed by a fixed number of the innermost 
stars from the initial disc. In the next step, we exclude from 
the CWS all the stars whose angular momenta deviate from 
the mean angular momentum of the CWS by more than 20° . 
On the other hand, the stars initially from outside the CWS, 
which do not fulfil the latter condition, are included into the 
CWS. Then, we recalculate the mean angular momentum of 
the CWS and repeat the whole procedure iterativelly until 
there are no changes of the CWS in between two subsequent 
steps. 

Observations sugg est (|Paumard et al.l 120061: 
iBartko et all 12009] , l2010l ) that roughly one half of the 
WR/OB stars are members of the observed CWS. We thus 



4 DISCUSSION 

As we have shown in the previous section, all the WR/OB 
stars may have been formed within a single gaseous disc 
and, subsequently, brought to their present location by the 
combined effects of two-body relaxation and differential pre- 
cession. In the following we establish a set of parameters 
for which the evolution of the young stellar system leads to 
a configuration compatible with the current observational 
data. For this purpose, we follow the evolution of the sys- 
tem for various values of the model parameters. Within the 
results, we then concentrate on Ncws/N, iews, (eews) and 
(e) and confront their values at t = 6 Myr with the obser- 
vations. 



4.1 System parameters compatible with the 
observations 

To begin with, we follow the evolution of the young stel- 
lar system for various values of A/cnd and M c with the 
other parameters fixed to their canonical values (see Table 
[TJ. The results indicate that the strongest constraints on 
the possible values of Mgnd and M c come from Ncws/N 
and iews- Their values for M c £ [0.01 M., 3 M,] are de- 
picted by the solid lines in the top panels of Fig. [3] whereas 
Mcnd remains constant along each of the lines and is set 
to either 0.1 M, or 0.3 M. or 0.6 M. . According to these re- 
sults, the evolution of the young stellar disc leads to values of 
Ncws/N which accommodate the observational constraints, 
if 0.1 M. < A'/cnd < 0.3 M. and 0.01 M. < M c < 2M.. 
However, the upper limit for the mass of the cluster of late- 
type stars has to be reduced to M c < M, since larger values 
do not lead to the observed iews ~ 90° . Both observational 
criteria are, therefore, fulfilled if 0.1 M, < M C nd < 0.3 M. 
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Figure 3. Number of stars within the CWS (top-left), inclination of the CWS with respect to the CND (top-right) and mean eccentricity 
of all the stars in the young stellar system (bottom-left) and within the CWS (bottom-right) at t = 6 Myr for various values of Mcnd 
and M c . All results are averaged over 12 realisations. Solid lines: ft = 7/4, N = 200, m G [4Mq, 120 Mg] (only properties of stars with 
m > 12 M Q displayed), a = 1. Dashed line: M C nd = 0.3 M,, /S = 1/2, N = 200, m G [4M , 120 Mq] (only properties of stars with 

/3 = 7/4, TV = 136, single mass, m = 50 Mg. Common parameters for all 



m ^ 12 Mq displayed), o = 1. Dotted line: 
models are: M. = 4 X 10 6 M Q , A = 2.5°, 



M CND = 0.3 Mm 
_ 7no 
CWS — ' u • 



and 0.01 M, < M c < M, . Moreover, the results of our com- 
putations with M c = show that even in this case, the 
evolution of the young stellar system leads to a configura- 
tion which matches the observational data (due to logaritmic 
scale in Fig. [3] the corresponding values are not displayed). 
Hence, we find the final intervals 0.1 M, < A/cnd < 0.3 M, 
and < M c < M.. If we substitute M, = 4 x 10 6 M©, the 
intervals transform to 4 x 10 5 M e < Mcnd < 1-2 x 10 6 M 
and < M c < 4 x 10 6 Mq . 

The intervals for allowed Mcnd and M c are not affected 
if we evaluate the eccentricity criterion. As can be seen in the 
bottom panels of Fig.[3] (solid lines), all the considered values 
of Mcnd and M c lead to similar values of both (eews) and 
(e), which satisfy the observational constraints. Neverthe- 
less, our results show that the orbital eccentricities in the 
young stellar system reach slightly higher values for lower 
M c . We attribute this effect to Kozai oscillations, which are 
less suppressed by the cluster of late-type stars. Further- 
more, our results suggest that also for M c ~ M. , all the 
orbits in the young stellar system gain somewhat larger ec- 
centricities, regardless the mass of the CND. Around the 



same value, lews appears to be more sensitive upon the 
variations of M c and Ncws/N reaches its minima (see the 
solid lines in the top panels of Fig. [3}. Hence, it seems that 
all of these effects are somehow connected with a stronger in- 
fluence of the CND on the dynamical evolution of the young 
stellar disc. However, at this point, we can not provide any 
explanation of this effect. 

In order to investigate whether the suggested inter- 
vals for Mcnd and M c depend upon the density pro- 
file of the cluster of late-type stars, we model the evo- 
lution of the young stellar system also for f3 — 1/2 and 
M c G [0.01 Mm , 3 Mm]- The other parameters remain at their 
canonical values (see Table [TJ. The dotted line in Fig. [3] 
proves that Ncws/N, as well as iews and both (eews) and 
(e), reach the same values as in the case with /3 = 7/4, 
except for the neighbourhood of the point M c ~ M. . The 
absence of the resonant effects observed in this case can be 
interpreted as a consequence of the different mass of the 
cluster of late-type stars enclosed within the young stellar 
disc, due to the different value of ft. However, since the value 
M c « Mm represents only the approximate upper boundary 
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of the suggested interval for M c , the latter effect is, for the 
purpose of this study, rather insignificant. 

Similarly, in order to test the dependence of the inter- 
vals upon the mass function of the young stellar disc itself, 
we perform a set of computations with the disc treated as a 
group of 136 single mass stars each with mass m = 50 M e . 
The other parameters are set to their canonical values. In 
this case, the total mass of the disc is the same as in all the 
models, which we have presented so far (ss 6.8 x 10 3 Mq). 
As demonstrated by the dashed line in Fig. [3] none of the 
results depend upon the mass function of the disc if its total 
mass is preserved. 

Our calculations further show that the evolution of the 
young stellar disc is not affected significantly if its total mass 
is changed within the range « 10 3 -10 4 Mq. On this account, 
all the results presented in the previous section would remain 
entirely unaffected even if we did not include the undetected 
stars with mass m £ [4 Mq, 12 Mq] as their overall mass 
represents only a small fraction of the total mass of the 
whole disc. 

On the other hand, decreasing the total mass of the disc 
by reducing the mass of the individiual stars m inhibits the 
combined effects of two-body relaxation and differential pre- 
cession. Namely, we do not observe the evolution of icws if 
m becomes negligibly small, i. e. if the stars in the disc can 
be considered as test particles. This effect is demonstrated 
by the bottom panel of Fig. [2 where we set m = 0.004 Mq. 
We can see that except for the few outermost stars, which 
are still slightly affected by Kozai oscillations caused by the 
CND, the inclination of the orbits remains constant through- 
out th e disc. Henc e , our results are in accord with the find- 
ings of ISubr et ail l|2009l) . 

We have determined the intervals for Mcnd and M c 
under assumption «cws — 70°. Our results indicate that 
if Jews ?J 60° and both the Mcnd and M c fall into the 
determined intervals, the evolution of the young stellar sys- 
tem leads within 6 Myr to a configuration that agrees with 
the current observations. With lower values of iews con- 
sidered, the CWS is entirely destroyed by the differential 
precession before it can reach the orientation perpendicular 
to the CND. On the other hand, considering iews closer to 
90° may increase the allowed intervals for Mcnd and M c . 
However, in order to fully understand the impact of different 
values of iews on t ne suggested intervals, a more detailed 
study would be required. We will focus on this issue in our 
future work. 

In this paper, we consider the young stellar disc to be 
initially thin. Our calculations show that the course of its 
evolution does not depend upon the value of its initial half- 
opening angle if Ao < 5°. 



4.2 Structure of the CND 

So far, we have modelled the CND as a single massive par- 
ticle on a circular orbit around the SMBH. This approxima- 
tion has been used instead of analytical descriptions, e. g. 
infinitesimally thin ring, for numerical reasons. Namely, due 
to its simplicity, single-particle approach minimizes the nec- 
essary computational time, and the corresponding perturb- 
ing particle can be implemented into the original NBODY6 
code in a trivial way. 

The single-particle approximation follows from the stan- 
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Figure 4. One of the realisations of the model with the CND 
treated as a group of A^nd = 20 equal-mass particles. The other 
parameters are set to their canonical values (see Table [TJ, except 
for Mcnd = 0.1 M. and M c = 0.1 M.. Top: Angular momenta 
of individual stars in the young stellar disc at t = 6 Myr (initial 
state denoted by an empty circle). Middle: Number of stars within 
the CWS (solid line). For comparison, we show the results for the 
canonical model (dashed line). Bottom: Eccentricity distribution 
at t = 6 Myr for all the stars in the young stellar system (empty 
boxes) and within the CWS (grey boxes). 
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dard averaging technique, which is commonly applied to 
many problems of celestial mechanics (see e. g. iMorbidellil 
2002). As a consequence, the single-particle approximation 
is equivalent to the model with the CND treated as an in- 
finitessimally thin ring if the assumptions of the averaging 
technique are satisfied. For the studied young stellar system, 
these assumptions can be written as the following two condi- 
tions for the orbital period P p of the massive CND particle: 
(i) P p must be significantly longer than the orbital periods 
of the early-type stars in the disc, and (ii) P p must be 
significantly shorter than the characteristic period P c of the 
studied phenomena. Since P c ~ 10 6 yr, P p ~ 10 5 yr, and 
P^ ~ 10 2 -10 4 yr, both conditions are fulfilled and, there- 
fore, the use of single-particle approximation is, in our case, 
well justified. 

The real CND is, rather than a ring-like structure, 
a gaseous torus, which consists of several somewhat au- 
tonomous clumps. In order to test whether the evolution 
of the young stellar disc can be affected by the dumpiness 
of the CND, we further consider the CND to be a group 
of iVcND equal-mass particles. It turns out that the CND 
constructed in this way is unstable with respect to its own 
gravity. Consequently, some of the particles successively mi- 
grate towards the SMBH. These particles can eventually be 
identified with several gaseous streams, which are indeed ob- 
served within th e radius of the CN D (for one of the most re- 
cent studies, see lZhao et al.ll2010l ). The results further show 
that all the combined effects of two-body relaxation and 
differential precession remain present (see Fig. [4] for the case 
with ./Vcnd = 20). Moreover, the infalling CND particles 
pose a stronger perturbation for the young stellar disc. As a 
result, the orbits in the disc gain higher eccentricities com- 
pared to models with the CND treated as a single massive 
particle on a stable orbit (see the bottom panel in Fig.[4] 
and the bottom-left panel in Fig.[l|. Similarly, the gradual 
deformation of the young stellar disc, as well as its eventual 
destruction, are also accelerated. 

Hence, it appears that the models with the CND treated 
as a single massive particle somewhat underestimate its in- 
fluence upon the dynamical evolution of the young stellar 
disc. On the other hand, the perturbative influence of the 
infalling parts of the gaseous CND would probably not be 
as strong as the impact of infalling point-like particles in the 
latter model. A more precise approach to the gas dynamics 
would thus be required in order to obtain a more accurate 
description of the CND. 



5 CONCLUSIONS 

We have modelled the orbital evolution of a thin self- 
gravitating stellar disc in a predefined external potential by 
means of N-body computations. In accord with the obser- 
vations of the Galactic Centre, we have considered the po- 
tential to include the gravity of the SMBH, the CND and 
the cluster of late-type stars. The results show that for a 
significant set of system parameters, the evolution of the 
disc leads to a configuration similar to that observed in the 
Sgr A* region within the lifetime of the WR/OB stars. In 
particular, while the outer parts of the disc are entirely dis- 
membered due to differential precession of the orbits caused 
by the CND, the inner parts remain undisturbed forming 



the CWS. Due to the influence of the CND, the CWS tends 
to change its orientation towards perpendicular with respect 
to the CND for a variety of initial configurations. Indeed, ac- 
cording to the observations, the CWS and the CND are in 
such a specific mutual orientation. It may thus be plausible 
that all the WR/OB stars observed in the Sgr A* region 
may have been born within a single gaseous disc and, subse- 
quently, brought to their present location by the combined 
effects of two-body relaxation and differential precession. 

We have further determined the possible values of the 
mass of the CND, Mcnd, which lead to a configuration 
compatible with the observational data. Based on several 
observational criteria, we have found the approximate in- 
terval: 0.1 M. < Mcnd < 0.3 M.. Although this interval 
represents a very strict constraint on the mass of the CND, 
it is in accord with the m ost recent observat i onal e stimate 
Mcnd « 10 6 M found bv lChristopher et al l d2005h. More- 
over, in agreement with the conclusions of ICuadra et al.l 
(2008), the suggested interval for Mcnd proves that unper- 
turbed two-body relaxation of the young stellar disc does 
not lead to the observed configuration of the WR/OB stars. 

Analogously, we have found the interval for allowed val- 
ues of the mass, M c , of the cluster of late- type stars enclosed 
within the radius of the CND: < M c < M.. Unlike the 
CND, this interval is rather wide and contains, from the ob- 
servational point of view, virtually all reasonable values of 
M c . We have confirmed that none of the suggested intervals 
depend neither upon the shape of the stellar mass function 
in the young stellar disc nor the density profile of the cluster 
of late-type stars. 

Let us note, however, that a less symmetric, e. g. flat- 
tened, cluster of late-type stars could have an impact on the 
evolution of the young stellar disc similarly to the influence 
of the CND. Nevertheless, since the current observational 
data do not provide any indication of such an asymetry, we 
have not considered this possibility in our analysis. More- 
over, due to the observed perpendicular mutual orientation 
of the CWS and the CND, it is likely that the CND indeed 
plays a crucial role in the evolution of the young stellar disc. 

Our results further indicate that the angular momenta 
of the majority of the WR/OB stars, which are scattered 
off the CWS plane, are expected to point to the same hemi- 
sphere with respect to the CND. This feature, together with 
the mass of the CND, as well as the mutual orientation of 
the CND and the CWS, represent strongly constrained out- 
put parameters of the scenario investigated in this paper. It 
would, therefore, be beneficial to concetrate observational 
efforts on the corresponding quantities and improve their 
accuracy in order to test the suggested hypothesis. 

Finally, let us emphasise that the combined effects of 
two-body relaxation and differential precession are efficient 
enough to transform the young stellar system into a config- 
uration which matches the observations in a period of time 
shorter than the estimated lifetime of the WR/OB stars. 
On this account, the presence of the CND in the system is 
not necessarilly required during the whole evolution of the 
young stellar disc. 
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